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Abstract 



Recently, effects of Feshbach resonances in atom-atom interactions were 
observed by varying the external magnetic field of an atomic Bose-Einstein 
condensate (BEC). We point out that the quasi-bound molecules created in 
the intermediate state of the resonance can form a second, molecular conden- 
sate. The many-body state of the system is then a hybrid atomic/molecular 
condensate with inter-condensate tunneling of atom pairs. A sudden varia- 
tion of the magnetic field results in oscillations of the number of atoms and 
molecules in their respective condensates, providing a signature of this novel 



1 



type of quantum tunneling. 
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Recently, Ketterle's group reported evidence for magnetic-field-induced Feshbach res- 
onances in the interparticle interactions of atomic-trap Bose-Einstein condensates 0- 
Near resonance, the two-body interactions are characterized by an effective interaction 
strength that can be tuned by varying the magnetic field, suggesting the construction of 
condensates with tunable interaction strength We present a theoretical interpretation of 
the condensate Feshbach resonances that predicts the formation of a molecular condensate 
that is coupled to the atomic condensate by quantum tunneling of atom pairs. 



I. THE BINARY COLLISION RESONANCE 

In an external magnetic field, the hyperfine induced spin flip of one atom in the presence 
of a second atom can bind the interacting atoms by bringing them to a spin configuration 
of higher energy. A resonance occurs if the potential of the bound channel supports an 



eigenstate that is nearly degenerate with the incident atoms 0- |]T0 . 

The detuning of the resonance, e, is the energy difference of the quasi-bound molecular 
state and the continuum of the incident atoms, pictured in Fig.(l). The energy difference 
A between the continua of the initial and quasi-bound channels, is the energy gain of a 
single spin flipped atom in the magnetic field B. When B takes on its resonant value, 
Bq, the energy of the quasi-bound molecule lines up with the continuum of the incident 
atoms so that A equals the binding energy of the bound state, A = Eb. Near resonance, 
A ^ Eb + (dA/dB) X [B — Bq], and one controls the detuning by varying the magnetic field, 
e = A-Eb^ (dA/dB) x [B - Bo]. 

In the many-body Hamiltonian, the spin flips to quasi-bound molecules are described by 

Hfr = a J d^r ipl^{r)'ipair)'ilJair) + h.c. , (1) 

where tprnji^m ( ^ajV'J) ^I's the annihilation and creation field operators of the molecules 
(atoms). The a-parameter in Eq.(|l|) is the transition matrix element proportional to the 
overlap of the molecular continuum and bound state wave functions. 
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The HpR interaction is responsible for the Feshbach resonance and gives a dispersive 
contribution to the effective interaction strength. To find the off-resonant detuning depen- 
dence we determine the energy shift of two atoms confined to a volume Q. If the initial 
state |ini) of energy Ei^i has both atoms in the same momentum state, then the interatomic 
interaction, [\a/2]J (Priljl{r)i(jl{r)-ijja{r)'ijja{r), contributes a shift AEa = A^/fi in first order 
perturbation. If we denote the quasi-bound state of energy Ei^t by |int), then the second or- 
der shift AEpR due to HpR is given by AEpR = | (int|ifFR|ini)p/[£'ini - Eint] = -[2a'^/e]/Q. 
Thus, AE = AEa + AEpR = Xes/^, where the effective interaction strength is equal to 
Aeflf = Aa - 2aVe- 



II. FESHBACH-RESONANCE INTERACTIONS IN THE CONDENSATE 

From momentum conservation, it follows that all quasi-bound molecules created from 
condensate atoms of vanishing momentum also occupy the same zero-momentum center-of- 
mass state, giving a first indication that the quasi-bound molecules form a condensate. A 
more insightful argument follows from the Heisenberg equations ihipa = [^Pa^H], ihtp^ = 
[ipmjH], where H is the Hamiltonian. The expectation value of these operator equations 
gives equations of motion for the condensate fields, 0^ = {ipm) and 0^ = {'4'a)- In the mean- 
field approximation (e.g. (■?/'a(r)^/'a(r)) ^ 0a(r)) fill, this procedure yields two coupled 
equations in place of the usual single condensate Gross-Pitaevskii equation [|l^- [Q: 

ifKpm = [ — + e + Xmrim + An^J^m + a4>l 

ih4>a = [ + KUa + \nm](t>a + 2a0*0rn , (2) 

where M denotes the mass of a single atom, n.^ and Ua represent the condensate densities, 
'^m = |0mP and Ha = l^aP, and Am, Aa and A represent the strength of the molecule-molecule, 
atom-atom and molecule-atom interactions. The a-terms that couple the equations describe 
tunneling of pairs of atoms between the (pm and 0a-fields. In particular, 0^ acquires a 
source term, a0^, so that its value cannot remain zero if (pa is finite: the tunneling creates 



a molecular condensate that is coherent with the atomic condensate. 



III. OFF-RESONANT STATICS 

Provided the condensate lifetime exceeds the time required by the system to reach its 
equilibrium, a discussion of the 'ground state' of the many-body system is meaningful. The 
equations that describe the static system are similar to Eqs.(2) with the time derivatives 
replaced by the chemical potential: itupa — ^ f^4>a and ih(j)fn '^f^4'm, where the chemical 
potential of the molecules is twice the chemical potential of the atoms, in accordance with 
the condition for chemical equilibrium. 

Of particular interest is the off-resonant (e > 0) limit: although e is tuned close enough 
to observe resonance effects, its value still exceeds the kinetic energy and interaction energies 
i^mnm, ^iT-m, Xa^a and An^). As a consequence, e >> fi, and << The ^^-equation 
yields (pm ~ —dUa/e (where we take (pa to be real) and the 0a-equation yields a Gross- 
Pitaevskii equation with the effective binary collision interaction strength: 



(Aa )na 



(3) 



2M ^ " e 

Note that equation by itself does not describe the appearance of the small molecular 
condensate, a feature that could have important applications |14 



In the off-resonant limit we can also understand the increase in loss rate that served as 
a signal to detect the Feshbach resonance. As the fraction of molecules remains low, we can 
assume that the atomic condensate decays mostly due to inelastic atom-atom collisions fia = 
—Caanl, and the molecular condensate due to molecule-atom collisions, = —Cmananm ~ 
—nlcmaio^/^]"^, where Caa and Cmm represent the corresponding rate coefficients. If we count 
each a-atom as one and each m-molecule as two condensate particles, the condensate particle 
density is n = Ua + 2n„^, and the condensate loss is described hj h = ha + 2hm ~ —n^{caa + 
2camn[a/ef). 

The observed Feshbach resonances create quasi-bound molecules of high vibrational 
quantum number. The single particle lifetime of such loosely bound molecules, limited 



due to collisions with atoms and other molecules, is substantially lower than the lifetime 
of the individual atoms. Estimates for the alkali-dimer rate coefficients give 10~^^ - 10~^ 
cm? I sec [|15|, compared to Caa ~ 10~^'*cm^sec~^. Thus, a purely molecular condensate of 
density lO^^cm"^ does not live longer than 10~^ seconds. In contrast, the molecular conden- 
sate in an off-resonant BEC survives much longer by compensating for the loss of molecules 
with atom pairs that tunnel in from the atomic condensate. 



The dynamics of a condensate in a time varying magnetic field, giving a time- dependent 
detuning e(t), is particularly interesting. A proper description requires an adjustment of 
Eqs.(2) to account for condensate loss. 

A formal treatment of the loss includes the channels of all chemical reactions that remove 
particles from the condensate. The two-body collision channels can be eliminated in per- 
turbation theory, modifying the equations of motion (Eqs.(|])) in a simple and predictable 
way: the interaction strengths become 'absorptive' with an imaginary part that determines 
the loss rates. In the 0rn-equation, for example, A ^ A — ifiCma/'^ and A^ — Am — ihcmmf^, 
where Cmm is the rate coefficient accounting for molecule-molecule collisions. 

In the off-resonant limit, the atomic density variations can be neglected and the atomic 
condensate serves as a reservoir of atoms. The appropriate approximations to the (pm- 
equation of Eqs.(|^), gives 



where 7^ is the molecule loss-rate: 7^ = Cmana- In the same off-resonant limit, 0^ ~ 

For the special case that the detuning undergoes a sudden shift to a value e/ and remains 
constant thereafter, we find 



IV. DYNAMICS 



m 
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where 0o is the initial field, 0o = 0m(^ = 0), and (p^o is the static off-resonant field value, 
0OO = — tt^a/[e/ + (A — 2Aa)na — ihjm/'^] ~ —dUa/ef. Note that the molecular density 
has an oscillating contribution ~ 2|(0o — 0oo)0^l cos([e/ + An^ — 2Aana)]t//i) exp(— 7„t/2). 
The molecules that appear and disappear during the oscillations tunnel in and out from 
the atomic condensate so that the atomic density has a contribution that oscillates out of 
phase with twice the amplitude. The oscillation is a genuine quantum tunneling effect and 
stems from the interference of the propagating initial field amplitude with the propagating 
amplitude of atom-pair tunneling. In general, the oscillations appear if the detuning was 
changed at a rate |e| that exceeds 167^1, or \B/[B — Bq]\ >> 7^- In the opposite limit, 
|e/e| < 7^, the system adiabatically follows its ground state. Experimentally, the oscillation 
can be observed, for example, by illuminating the BEC with light that is resonant with a 
transition of the quasi-bound molecule. The intensity of the image will be modulated at the 
frequency of the quantum tunneling oscillations. 

Near resonance, the nonlinear density variations cannot be neglected, but the dynamical 
behavior remains qualitatively similar: a sudden shift of detuning results in out-of-phase 
oscillations of the atomic and molecular condensate densities, as shown in Fig. (2). The 
oscillations are damped and the condensate density decays on a longer time scale. 



In an atomic BEC, the molecules created in the intermediate state of a Feshbach reso- 
nance form a second condensate that is coherent with the atomic condensate field. Pairs of 
atoms tunnel between the condensates. We have discussed a clear signature of this novel type 
of quantum tunneling: a sudden change in the external magnetic field results in oscillations 
of the number of atoms and molecules in their respective condensates. 
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Figure Captions 

Fig.l : Schematic representation of the molecular potentials of the incident and intermediate 
state channels . The energy difference of the continuum levels, A, is the sum of the binding 
energy Ei, of the quasi-bound state and the 'detuning' e. 

Fig. 2 : Plot of the particle densities: the total condensate density, n — na-\- ^rim, in full 
line, the atomic density in dashed line and the molecular density in dash-dotted line. 
The oscillatory behavior of the atomic and molecular densities is a signature of quantum 
tunneling. The calculation is for a homogeneous BEC that was initially in equilibrium at 
density n — lO^^cm"^ when the detuning experienced a sudden shift from e = 50An to 
e = 2Xn. The order of magnitude of the interaction parameters, An = A^n = A^n = 
ay/2n — 10^ Hz, and of the decay parameters, c^a = Cmm — 5 x 10~^^cm^sec~^ (while 
neglecting the atomic decay) are reahstic. 
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